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Abstract 

Starting with the quasi-Bell states of the qubit-oscillator system, we obtain time evolution 
of the density matrix under the adiabatic approximation. The composite density matrix 
leads to, via partial tracing of the qubit degree of freedom, the reduced density matrix of 
the oscillator that is utilized to obtain the quasi-probability distributions such as Glauber- 
Sudarshan P function, Wigner W function and Husimi Q function. The negativity of the 
Wigner function acts as a measure of the nonclassicality of the state. The negativity becomes 
particularly relevant in understanding a comparison between the Wigner entropy with the 
Wehrl entropy, which are based on the W function and Q function, respectively. 


I Introduction 

A two-level system (qubit) that interacts with a radiation field represented by a single oscil¬ 
lator mode is one of the important models in quantnm optics. This model has been studied 
extensively under rotating wave approximation (RWA) [T] that holds good for weak coupling 
regime. Following the recent developments in experiments [2][3], it is of interest to study the 
system with higher qubit-oscillator coupling. In the (ultra) strong coupling domain and for low 
qubit frequency where the oscillator frequency sets the scale, the system has been investigated 
using adiabatic approximation that utilizes the separation of the time scales of the high 

oscillator frequency and the low (renormalized) frequency of the qubit. 

To study the dynamics of the quantum state, we employ various quasi-probability distri¬ 
butions on the phase space as preferred tools. For instance, the Wigner function [7] gives a 
connection between the classical and qnantum dynamics. Unlike the classical (true) probability 
distributions, the Wigner function W may assume negative values. This negativity [H] serves as 
an indicator of the nonclassicality of the quantum state. On the other hand, entropy is a key 
concept for the quantum systems that provides a framework for studying the loss of information. 
The W distribution has also been recently employed mm for developing a new measure of 
quantum entropy. This is found to have close correspondence with the negativity parameter. 

The present manuscript is organized as follows; The initial quasi-Bell composite state and 
the corresponding reduced density matrix of the oscillator are considered under adiabatic ap¬ 
proximation in Sec.II. This density matrix is utilized in Sec.Ill to derive the Glauber-Sudarshan 
P distribution mm- In Sec.IV we evaluate the Wigner W fnnction [7] which is used to study 
the nonclassicality and decoherence of the system. The relevant Husimi Q distribution [S] ad¬ 
mits, in the weak coupling regime, a closed from evaluation d la the procedure employed in 
[S]. Comparison of the time evolution characteristics of the Wherl entropy [H] based on the Q 
distribution, and that of the Wigner entropy mm, which is based on the W distribution, is 
discussed in the context of negativity of the associated quantum state. 


I 


II Reduced density matrix of the oscillator 


We study a coupled qubit-oscillator system with the Hamiltonian that reads in natural 

units (h = 1) as follows: 



-I- a + Adz (a^ + a), 


( 2 . 1 ) 


where the harmonic oscillator with a frequency a; is described by the raising and lowering 
operators {a\ a\h = a^a). The qubit characterized by an energy splitting A as well as an 
external static bias e is expressed via the spin variables (dx, dz)- The qubit-oscillator coupling 
strength is denoted by A. The Fock states {h|n) = n|n), n = 0,1,... ; a |n) = ^/n |re—1), |n) = 

\/n + 1 \ n + 1)} provide the basis for the oscillator, whereas the eigenstates + 1) = + | + 1) 
span the space of the qubit. The Hamiltonian ()2.ip is not known to be exactly solvable. In 
the present work we follow the adiabatic approximation mm that hinges on the separation 
of the time scales governed by the high oscillator frequency and the (renormalized) low qubit 
frequency: w » A. 

We consider the following quasi-Bell initial states of the coupled system 

|'!/;(0))^-^ = —p (|1, a) + I — 1, —a)), |a) = D(q;)| 0), D(a) = exp(Q;a^ — a*a), (2.2) 

v2 


where the coherent state {|a) Va = Re(a) -I- Hm(Q;) e C} for the oscillator degree of freedom is 
realized by the action of the displacement operator on the vacuum. Under adiabatic approxi¬ 
mation [5], the evolution of the initial state (12.211 yields the bipartite density matrix: 

pl«(i) - (2.3) 


The partial tracing, say, over the qubit degree of freedom produces the reduced density matrix 
of the oscillator: 

(2.4) 

that may be expressed US], in terms of the displaced number states |n+) = (+A/a;) |n) as 

follows: 


1 f \n f ^ \TTi 

= 9exp(-|a+|2) 2 r-r^ — Uiit) ^mity\n+Xm+\ 

2 Vnlml V 


n,m=0 


+(-l)”’^™C(^)* exp ( - i(n - m)ut). (2.5) 


The coefficients in (j2.5p are given by 

^n{t) = exp(ix„t) -t exp(-ixnt), ^ ^ ^ ^ X X ^ 

2Xn 


Xn - e + 
2Xn 


where, Xn = V'^n + X = -—Ln{x), x = (2A/w)^ , A = Aexp (-§) : « = | and a+ = 
a + X/lo. The reduced density matrix (12.5p of the oscillator may be utilized to calculate the 
quasi-probability distributions. 


Ill Phase space distributions 

A. Glauber-Sudarshan P function 

The well-known Glauber-Sudarshan P function [13 [S] admits a diagonal representation of the 
oscillator density matrix in the coherent state basis: 

p = ^pmw\d^x (3.6) 
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For an arbitrary quantum state the relation (I3.6p may be inverted and the P function is uniquely 
expressed as 


^ (3.7) 

For our choice of the reduced density matrix (I2.5p the P function may be derived readily as 
follows: 


F(±>(/3) - -exp(-|a+|=) ^ 


n,m=0 


X I C4: I 

nlml 




Vd/ 3 / ' 

<d/ 3 */ 




where /3+ = /3 + ^. The above P distribution incorporating derivatives of <5 functions is highly 
singular. This is a typical behavior observed for nonclassical states. Due to its highly nonsin¬ 
gular nature, employing P{f3) distribution directly towards producing a quantitative measure of 
nonclassicality is complicated. So other quasi-probability distributions should be considered in 
this regard. 


B. The Wigner W function 

For an arbitrary density matrix p the Wigner quasi-probability distribution is defined [7] via the 
displacement operator as 

VF(/?) = J Tr[pD(7)] exp(/37* — /3*7)d^7. (3.9) 

But the evaluation of the Wigner function using the definition (13.9p is not always easy. An 
alternate series representation of the distribution W{j5) in terms of the diagonal matrix elements 
in the displaced number states is known |16j : 

2 ^ 

W{/3) = - V (-l)\/3, k\p\/3, k), 1/3, k) = D(/3)|A:>. (3.10) 

^ to 


Substituting (|2.5p in (|3.1Up we obtain the time evolution of the Wigner function for the initial 
quasi-Bell states: 


1 


. -exp(-|a+7X;{-l)‘ 

vr rt. 


k=0 


2j -Tw- Dmk{/3+)e 

m=0 


iuimt 


77 irv* i™(r+iti 


m=0 


\/mJ. 


(3.11) 


The component DmA:(/3) is evaluated using the hypergeometric function as 




D^fc(/3) = (-!)" exp(-|/3|72) 


VudM 


Fo -m,-k;_;- 


l/3| 


2 / ’ 


(3.12) 


where the generalized hyper geometric function 2 Tb(m, n;_; z) reads 


CC / 

2Foim,n]_]z) = ^ - 


m)r{n) 




r=0 


t\ 


(3.13) 
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Figure 1: The density plots of using (I3.15P for A = 0.15ti;,e = 0.03a; and a = 3 at 

various values of scaled time ujt (a) 0, (b) 50, (c) 100 for the coupling strength A = 0.3 a;. 


The identity 

/ 1 \^/, 1 \ 

- 2^0 \y-k, -m; -j 2 Tb (^ - fc, -n; -j 

k=0 

yields the Wigner function (|3.11l) in the following form: 



= 2 ^ 


,n+m X 


e '' 2 -Fo(^ - m, -n;_; 



(3.14) 


VF(±)(/3) 


1 * (‘In i™ 

= -exp(-|a+p) ^ —exp(-ia;(n - m)t) [exp(-2|/3+|^) x 

TT ^ r. n\m\ ^ 

n,m=0 

X (/3:)"(/3+)™e:±(t)C±(t)^Fo + (-1)"+™ x 

X exp(-2|^_|2)(^*)"(^_)™(J:+(t)*e:+(t)2Fo . 


(3.15) 


The expression ()3.15l) satisfies the normalization condition: 


(3.16) 

FigH] shows the Wigner function at three different instances of time. The distribution for the 
initial state consists of two Gaussian peaks as it is the superposition of two well-separated 
coherent states. As time evolves it assumes negative values demonstrating the nonclassical 
nature of the state. For a suitably strong qubit-oscillator coupling a large number of interacting 
modes set in. The quantum interference between these modes give rise to negative values of the 
VF-distribution in the zone of the phase space intermediate between the positive peaks. The 
volume of the negative domain of the Wigner function on the phase space is considered as a 
quantitative measure of nonclassicality of the density matrix [8]: 

A(±) = J|w(±)(,5)|d2^-1. (3.17) 

A non-zero value of A^-) indicates the characteristic quantum nature of the given state. The 
time evolution of A^-^ for various coupling regime is discussed in the following subsection. 
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C. Husimi Q function 

The Husimi Q function [9] is defined as the expectation value of the density matrix of the 
oscillator in an arbitrary coherent state: 

Qi/3) = -</?|p|/3>. (3.18) 

TT 

For the state f)2.5p under consideration the Q-function assumes the explicit positive definite form 


Q^-\f3) = —exp(-|a+|2)(exp(-|/3+|2)|X±|2+exp(-|/3_|2)|2}± 
The Fourier sums in (13.191) are given by 


(3.19) 


+ ^ (a+/3*)"' ^ /■ ■ to nr>\ 

T-= -j-e:^(t) exp(-^na;^), (-1) - exp{-mu}t). (3.20) 

n=0 n=0 


n\ 


The Husimi distribution (j3.19p of the reduced density matrix (j2.5p does not have any zero on 
the phase space except at asymptotically large radial distances. It satisfies the normalization 
condition: (3 = 1. Adopting the procedure developed in [5] where the Laguerre 

functions are approximated by their linear parts in the regime ^ < 0.1 we may approximate the 
Fourier sums (13.201) in closed forms as 


^(±) = 


2 ) 


>1) cos (^$ 4 ^ - 0 + iexp(-^>|) sin 

exp($!) sin + x sin j 

,2 . , , ... 


exp(<I>] 
, e 


(+) = 


exp(-$|) sin + ^) “ sin 

0 + *exp(Tl) sin ^ 


+i 

exp(—'I'l) cos 

. e 

—i— 

A 


T 


exp(— tJ) sin ~ ('^t^ 


(^ ■ 


(3.21) 


In ()3.21l) we have used the notations = a+/3+exp(—fcjt), = a+/3* exp(—fcjf), = 

XT II XT XT XT ^ 

cos —, TJ = 'I't cos —, = <I>t sin —, sin — and t = At. 

Employing the following well-known interrelations m the Wigner and Husimi functions can 
also be directly obtained from the P distribution: 

g(/3) = 1 J e-l^-T'lV(7)(i2^, W{P) = i J e-2|^-T'lV(7)(i2^. (3.22) 

On the other hand, the following property HZj suggests that the Q function may be considered 
as a ‘coarse-grained’ behavior of the IF function: 

W) = ^Je"''^"^'V(7)A, (3.23) 

where the Q function is obtained after a suitable ‘smearing’ of the IF function with a positive 
definite kernel. The results obtained earlier for the quasi-probability distributions using the 
direct formulae and the interrelation are found to be consistent. To understand the role of 
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Figure 2: The time evolutions of negativity (red), Wigner entropy 5^^ (green) and Wehrl 
entropy obtained for the values of A = 0.15w,e = 0.03a; at (a) a = 2, A=0.08a; and (b) 
a = 3, A= 0.3a;. 


negativity of the Wigner distribution we now do a comparative study of the Wehrl entropy [2] 
that is based on the Q function and the recently proposed m measure of entropy based on the 
Wigner W function. The Wehrl entropy [T^ defined as 

= - JlogQ(±)(/3)d2/3 (3.24) 

acts as an information-theoretic measnre describing the delocalization of the oscillator on the 
phase space. It is also of interest to study the quantum entropy based on the modulus of the 
Wigner distribution |IT(/3)| |11] which is a non-negative quantity: 

^ = - J \W^-\P )\log |IT(±)(/3)|d2/3. (3.25) 


It is evident from the Figs. [2] (a) and (b) that the time evolutions of the Wigner entropy 
(I3.25P and negativity parameter (|3.17p have close kinship with each other. An increasing 5^^ 
reciprocates increment in A^'^^ and vice versa. In this sense the Wigner entropy 5^^ reveals the 
extent of nonclassicality of a quantum density matrix. We distinguish between two possible sce¬ 
narios depending upon the qubit-oscillator coupling strength, (i) In the strong coupling regime 
(A/o; < 0.1) we observe (Fig. [2] (a)) that a periodic structure in the phase space distribution 
develop, and, conseqnently, the physical variables such as the Wigner entropy -5^^ negativity 
and Wehrl entropy follow similar periodic patterns that may be identified with the 
revival and collapse of the qubit density matrix elements. For a dominant value of when 
quantum interference effects are overwhelming, we, expectedly, find On the other 

hand, for a low A*-^) region the inequality is reversed: < *5'^^. As the Q-function is ob¬ 

tained from the IT-fnnction (I3.22p after suitable averaging with a positive definite kernel, the 
Wehrl entropy <5^^ would be larger in a low negativity A^'*'^ regime [10] . (ii) In the ultra strong 
coupling regime (A/o; » 0.1) all Fourier modes for the qubit-oscillator interaction are excited and 
a fully randomized interference pattern very quickly evolves. The randomized interferences of 
a large number of Fourier modes lead to quasi-stationary values of the phase space observables 
(Fig. [2] (b)). Moreover, these interferences necessarily develop areas on the phase space with 
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negative values of Vl^-distribution. The average value of the negativity parameter increases 
with increase in the coupling strength. To keep, however, the normalization sum rule (|3.16p 
intact, suitable increment in the magnitude of the W function takes place leading to increased 
value of the entropy Sw- In the quasi-stationary state the negativity is statistically pre¬ 
served. The quasi-stabilization of A^^) occurs after a suitable decoherence time. Consequently, 
we observe that except for a brief initial period the Wigner entropy is consistently more 
than the Wehrl entropy even though the positive definite Q function may viewed as the 

smeared form of the W distribution. 

IV Conclusion 

Phase space dynamics of the strongly coupled qubit-oscillator system is studied under adiabatic 
approximation. Using this approximation the oscillator density matrix have been written in 
terms of the displaced number states. This density matrix is employed to calculate the quasi¬ 
probability distributions such as Glauber-Sudarshan F function, Wigner W function and the 
Husimi Q function. The negative values of the Wigner distribution acts as an witness of the 
nonclassicality of the state. The time evolution of this nonclassicality parameter of the state is 
obtained. In the ultra strong coupling regime the negativity assumes, after a suitable decoher¬ 
ence time, a quasi-stationary value. The quantum entropy based on Wigner function exhibits 
qualitatively similar behavior as that of the nonclassicality. The value on the nonclassicality 
measure serves as the key to understand the relative values of the entropies based on Wigner 
and Husimi distributions. 
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